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Abstract. We consider a type of intermittent behavior that occurs as the result of the interplay between 
dynamical mechanisms giving rise to type-II intermittency and random dynamics. We analytically deduce 
' . the law for the distribution of the laminar phases, which has never been obtained hitherto. The already 

I known dependence of the mean length of the laminar phases on the criticality parameter [PRE 68 (2003) 

. ^ • 036203] follows as a corollary of the carried out research. We also prove that this dependence obtained 

earlier under the assumption of the fixed form of the reinjection probability does not depend on the 
relaminarization properties, and, correspondingly, the obtained expression of the mean length of the laminar 

> 

00 ' phases on the criticality parameter remains correct for different types of the reinjection probability. 

o 

' PACS. 05.45.-a Nonlinear dynamics and nonlinear dynamical systems - 05.40.-a Fluctuation phenomena, 

o ■ 

■ random processes, noise, and Brownian motion 



■ 

X 



1 Introduction nonlinear systems. In particular, it was discovered that 



^ ' random fluctuations can actually induce some degree of 

Intermittency is known to be an ubiquitous phenomenon, i ■ i ■ ^ r i- ^ irnm-nrm i 

-1 : order m a large variety or nonlmear systems |91 I10II11) , and 

with its arousal and main statistical properties having , , -111 i • 1 ^ i • 1 

° sucn phenomena were widely observed m relevant physical 

been studied and characterized already since long time . , fT^n-rn 

" circumstances |61I12). 

ago. The different types of intermittency have been classi- 
fied as types I-III [I], on-off intermittency eyelet There are no doubts that different types of intermit- 
intermittency [5J|5] and ring intermittency [71IH] . From the tent behavior may take place in the presence of noise and 
other side, increasing interest has been put recently in the fluctuations in a wide spectrum of systems, including cases 
study of the constructive role of noise and fluctuations in of practical interest for applications in physical, radio en- 
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gineering and other applied sciences. It is plausible that 
such an interaction would originate new types of dynam- 
ics. Therefore, the intermittent behavior in the presence 
of noise has been studied by means of Fokker-Plank equa- 
tion |13j and adopting renormalization group analysis [14 , 
but the characteristic relations were obtained only in the 
subcritical region, where the intermittent behavior is ob- 
served both in the presence of noise and without noise. 
Recently [TCllTFlfTT] , the theoretical and experimental con- 
sideration of the intermittent behavior in the presence 
of noise has been considered in the supercritical region 
(where intermittency does not take place in the absence 
of noise) for the type-I intermittency. 

Obviously, in the presence of noise the other types of 
intermittent behavior mentioned above may result in the 
distinct types of dynamics. In this paper we report for 
the first time the important characteristic (namely, the 
distribution of the laminar phase lengths) deduced ana- 
lytically for the type-II intermittency in the presence of 
noise for the region of the supercritical parameter values, 
which has never been obtained hitherto. This characteris- 
tic is of great importance, since scientists studying the in- 
termittent behavior of the complex systems often can not 
variate the criticality parameter in the experiments (e.g., 
in physiological and biological systems), and, therefore, 
they can deal only with the the distribution of the lam- 
inar phase lengths, whereas the dependence of the mean 
length of the laminar phases on the criticality parameter 
can not be obtained. The already known dependence of 
the mean length of the laminar phases on the criticality 
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parameter 18 follows as a corollary of the carried out re- 
search. Moreover, we prove that this dependence obtained 
in |18| under the assumption of the fixed reinjection prob- 
ability (taken in the forms of delta-function and uniform 
distribution) does not depend practically on the relami- 
narization properties, and, correspondingly, the obtained 
expression of the mean length of the laminar phases on the 
criticality parameter remains correct for different forms of 
the reinjection probability. The obtained analytical distri- 
bution of the laminar phase length is verified by means of 
numerical calculations of the model system dynamics. 

2 Analytical approach 

The standard model that is used to study the type-II in- 
termittency |18) is the one-parameter cubic map 

Xn+l = .f{Xn) = {I + e)Xn + xl, (1) 

where e is a control parameter. Below the critical param- 
eter value (i.e., for e < ec), the stable fixed point Xc — 
is observed, while above Ec this fixed point Xc = be- 
comes unstable and the point representing the state of 
the map ^ moves around it with slowly increasing am- 
plitude. This movement in the vicinity of the fixed point 
corresponds to the laminar phase, its mean length T being 
inversely proportional to (e — Ec), i.e. 

r^(6-6,)-i. (2) 

To develop the theory of type-II intermittency in the 
presence of noise, we consider the same cubic map ^ with 
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the addition of a stochastic term ^„ 

Xn+l = f{Xn) = (1 + e)Xn + + (3) 

where ^„ is supposed to be a delta-correlated white noise 

[(6.) =0, i^nU) ^DSin-m)]. 

The influence of the stochastic term ^„ on the behav- 
ior of the system is governed by the value of parameter D. 
For positive values of the control parameter e (e > 0), the 
point corresponding to the behavior of system ([3]) moves 
in the iteration diagram around the unstable fixed point, 
its motion being perturbed by the stochastic force. As far 
as the intensity of the noise is not large, the characteris- 
tics being close to the classical type-II intermittency are 
observed. 

A different scenario occurs for control parameters e as- 
suming negative values (e — —e, where e = |e| > 0). In this 
case, the point corresponding to the behavior of system ([3]) 
is localized for a long time in the region — ^/s < x < \fe 
and its dynamics is also perturbed by the stochastic force. 
As soon as the system state point arrives at one of the 
boundaries Xh — ±£^^/^ due to the influence of noise, a 
turbulent phase arises, though such kind of events is very 
rare. 

In this case, the behavior of the map ^ differs radi- 
cally from the dynamics of the system ((Ij, since the tur- 
bulent phases are not observed for e < if there is no 
noise. Therefore, such a region of negative values of the 
e-parameter is the main subject of interest for the type-II 
intermittency in the presence of noise. 
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2.1 Probability density 

Having supposed that: (i) the value of e is negative and 
rather small and (ii) the value of x changes per one it- 
eration insufficiently, we can consider — Xn) as the 
time derivative x and undergo from the system with dis- 
crete time to the flow system, in the same way as in 
the case of the classical theory of the intermittency [T] . 

Since the stochastic term is present in (jS]) we have to 
examine the stochastic differential equation 

dX = X{X'^ - e) dt + dW (4) 

(where X(t) is a stochastic process, W(t) is a one-dimensional 
Winner process, e = |e|) instead of the ordinary differen- 
tial equation dxjdt = xix? + e) considered in the classical 
theory of type II intermittency. 

The stochastic differential equation (|H) is equivalent to 
the Fokker-Plank equation 

for the probability density px{x,t) of the stochastic pro- 
cess X(t). The chosen initial condition is px{x, 0) = S{x — A), 
where 6{-) is a delta- function. Such a choice of the initial 
form of the probability density px{x,0) corresponds to 
the beginning of the laminar phase, when the point rep- 
resenting the state of the system ([3]) is in the place with 
coordinate x — A {—e^^^ < A < e^l'^) at time t = 0. In 
other words, we suppose that the reinjection probability 
is a (5-function 

P,,,(x) = 8{x - A) (6) 

and after the relaminarization process the system is al- 
ways returned to the state x ~ A. Although the rein- 
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jection probability Pin{x) is well-known to be important As the maximum of the probability density should coin- 
factor and should be taken into account when the statisti- cide with the stable fixed point = 0, one has 
cal properties of the intermittent behavior are studied 1191 

9'{z) = 0. (9) 

120) . in the considered problem the form of the reinjection 

probability practically does not influence on the distribu- Under the mentioned assumption, we consider the or- 
tion of the laminar phase lengths (and the dependence of ^^^^^^ differential equation 

the mean laminar phase length on the criticality parame- D* g"{z) — 2 (z(z^ — l)g(z))' — (10) 

ter, respectively), as it will be shown below. 

instead of Q for the region z e (— 1; 1). 

To reduce the number of the control parameters the 

This equation is equivalent to 

normalization z — t — te may be used, after which 



Eq. ® may be rewritten in the form D*g'{z) - 2z{z - l)g{z) + Ci - 0, (11) 

where Ci is constant. To solve this equation we use the 
Ot dz"- '2 dz^ integrating factor 



z2 



where D* = De-^ pziz,T) = pxize^/^TE-^). ^^^^ ^ ''''^ [d^ ~ Y ) ) ' ^^^^ 

As the coordinate of the system state stays for a long xhe solution of ^ may be found in the form 
time in the region |z| < Zf, = 1, one can suppose that the 



Jo ''''^[d^V'Y)) '^'^^^ 
probability density may found the form of the metastable 5(^) — — 7 ^2 \ \ ■ ^^^) 

D* exp 



distribution decaying slowly for a long period of time. \D* 

n r ^1 1 J i • From Eq. ([T31) one can obtain easily that g'(0) = 2C1/Z?*. 

Ihe relaxation process of the probability density to this ^ — ^ ■' ^ \ > '-1 

^,,1^^. Ti- Taking into account the condition (191), one comes to the 

metastable state is supposed to be very fast m compar- ^ 

conclusion that Ci = 0. Note, in this case the obtained 



function 



ison with the time of the metastable distribution decay, 
therefore, one can neglect the transient Q <t <ttr- Under 
the assumptions made above the probability density may ^( _ 

be written in the form pz{z, r) = A{T)g{z)^ Vz e (^1; 1), 
where A{t) > decreases very slowly as time increases, 
i.e. dA/dr w 0. The function g{z) should satisfy the con- 
ditions pz{z, t) ~ A{t) exp ( --^ ( 1 - y ) ) • (15) 



g(z) = -^exp -— 1-y ) ). (14) 



also satisfies the conditions ([5]). Therefore, the probability 
density pz{z,t) in the region z G (—1; 1) is 



1 

g(z) > Vz e (-1; 1) and J g{z)dz< 00. (8) 

-1 



When the laminar phase is interrupted, the system es- 
capes from a metastable state. Therefore, we suppose that 
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the decrease of A{t) should be determined by the proba- 
bility distribution taken in the boundary points z — ±1, 
i.e., dA{T)/dT ~ —pz{±l,T). This assumption, which is 
also equivalent to neglecting the time correlation of the 
orbit, may be rewritten as 



dAjr) 
dT 



-fcy4(r) exp 



1 



2D* 



(16) 



where is a proportionality coefhcient. Evidently, the de- 
crease of A{t) is described by the exponential law 

A(t) = yl(0)exp(-fc7?T), r] = exp{~l/{2D*)). (17) 

Having returned to the initial variables x and t we 
derive the following expression for the probability density 



px{x,t) ~ A(t)exp ( -— ( e - Y 



where 



and 



A(t)=A(0)exp( -- 



(18) 



(19) 



(20) 



with A{t) being considered as a normalizing factor, i.e., 
A{0) J exp(^-^(^e-^)) dx ^ 1. (21) 

To confirm the assumptions made above and the ob- 
tained equations, we have compared the evolution of the 
probability density px{x,t) given by (jlSp with the re- 
sult of the direct numerical calculation of the Fokker- 
Plank equation ([5]) with the values of control parameters 
e = 10"^ D = 10"^ 

The evolution of the probability density px{x,t) ob- 
tained by the numerical calculation of ([S]) is shown in 




500 



a 




b 

Fig. 1. The evolution of the probability density px(x,t) 
obtained by means of the direct numerical integration of 
Fokker-Plank equation ©, e = 10"^, D = 10"^. (a) The ini- 
tial fragment of the density evolution involving the transient 
(0 < t < ttr, ttr ~ S X 10^). (b) The long-time evolution of 
px{x, t), with the transient being omitted, 5 x 10"^ <t<5x 10* 

Fig. [TJ One can see that after the very short transient 
< t < ttr the probability density px{x, t) arrives the state 
being close to stationary (Fig. [T]a). After that the value 
of px{x,t) decreases very slowly (according to the expo- 
nential law) with time increasing, with the form of the 
dependence of the probability density on x-coordinate be- 
ing invariable (Fig. [T]6). 

Fig. [5] also shows the profiles of the probability density 
px{x,t*) taken in the different moments of time. It is ev- 
ident, that after a very short transient (curve 1, = 15), 
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the density px{x, t) practically does not change when time 
increases. 

Two different profiles px{x^t*) corresponding to the 
time moments = 5 lO'^ and = 3 x 10'^ (cmves 2 
and 3, respectively) are very close to each other despite of 
the large time interval Z\t = tg — between them. More- 
over, they are in very good agreement with the approxi- 
mated solution A{Q)g{x) described by Eq. (fTS]) and shown 
in Fig. [2]by means of squares. As time goes on, the ampli- 
tude of the probability density decreases according to the 
exponential law, but very slowly (see Fig. [21 curves 4 and 
5, t\ — 10'' and = 5 x 10*, respectively), although the 
probability density form remains the same for all times. 

Therefore, taking into account the results of the di- 
rect numerical calculations of Fokker-Plank equation ([SJ 
and the comparison with the obtained approximated solu- 
tion (jl8p , we come to the conclusion that our assumptions 
are correct and can be used for the further analysis. 

The evolution of the probability density px{x,t) may 
be considered separately on two time intervals < t < ttr 
and ttr < t < +00, respectively. The first time interval cor- 
responds to the transient when the probability density 
px{x, t) evolves to the form (fT5)) being close to stationary. 
Only when < t < ttr the form of the reinjection proba- 
bility Pin{x) may infiuence on the evolution of the prob- 
ability density px{x,t). For t > ttr (when the transient is 
elapsed), the evolution of the probability density is defined 
completely by Eq. and it does not depend entirely on 
the reinjection probability Pj„(x). Since the transient is 
very short in comparison with the exponential decrease of 




Fig. 2. The profiles of the probability density px{x,t*) taken 
in the different moments of time t* obtained from the direct nu- 
merical calculation of Fokker-Plank equation ((S)). In the frame 
the dependence of the probability density p(0, t) in the station- 
ary point X = Q on time t (+) and its exponential approxima- 
tion (solid line) are shown in the logarithmic scale 

the probability density px{x,t) we can neglect them and 
use only the second time interval ttr < t < +oo to obtain 
the statistical characteristics of the type-II intermittent 
behavior in the presence of noise. It is clear, that in this 
case the obtained results do not depend on the rclaminar- 
ization process and the reinjection probability Pin{x). 

2.2 Distribution of the laminar phase lenghts 

Let us discuss now the relationship between the probabil- 
ity density px {x, t) and the distribution pit) of the laminar 
phase lenghts t. For this purpose we consider the ensem- 
ble of systems Let there be N{t) systems which at the 
moment of time t are in the stage of the laminar phase. Ob- 
vioulsy, the states x of these systems should be distributed 
over interval (— -v/e, \/e) according to the probability den- 
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sity px{x,t). After the infinitely small time interval dt 
several systems stop demonstrating the laminar behavior 
(the turbulent phase begins in these systems) and, cor- 
respondingly, their states x leave the interval (— -y/e, y^), 
while the states of other N{t + dt) systems demonstrating 
the laminar dynamics as before are distributed in accor- 
dance with the probability density px{x,t + dt)). If the 
total number of the systems in the ensemble under con- 
sideration is iVo, the number of the systems in which the 
laminar phase is finished during time interval [t, t + dt) is 

N{t) - N{t + dt) = 

f (22) 
= Nq [px{x,t) - px{x,t + dt)]dx. 

Since the number of the systems in the ensemble in which 
the laminar phase is interrupted at the moment of time 
falling in the time interval between t and t+dt is connected 
with the distribution p{t) of the laminar phase lengths t 
as 

N{t) - N{t + dt) = Nop{t) dt, (23) 
one can obtain the relationship 

Pit)^~j'-^dx (24) 

between the probability density px{x,t) and the distribu- 
tion p{t) of the laminar phase lenghts. 

Using relations (ITSl) , (fTO|) and (PT|) one can obtain, that 
the laminar phase distribution is governed by the exponen- 
tial law 

p{t)=T-'cxpi-t/T), (25) 

where T defined by Eq. ((20|) is the mean length of the lam- 
inar phases. The obtained expression (|20p for the mean 
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length T of the laminar phases is consistent with the for- 
mal solution derived in the previous considerations |18[ 
I21| . that may be considered as an additional evidence of 
the correctness of our results. Nevertheless, we would like 
to emphasize that the applicability of the results obtained 
earlier analytically [T51[^ are limited greatly by the as- 
sumptions (made to deduce corresponding equation) con- 
cerning the character of the reinjection process. Indeed, 
it is well known that the characteristics of the intermit- 
tent behavior even in the case without noise depend not 
only on the structure of the local Poincare map but on 
the reinjection probability distribution, with this depen- 
dence being sufficient (see e.g., [TOll^ ). At the same time, 
the analytical results mentioned above were obtained only 
under assumptions of uniform and delta-function reinjec- 
tion probability distribution. Our study allows to extend 
known relation to all types of the reinjection probability 
distribution, since based on the consideration carried out 
above we state that Eqs. (j20p and ([25]) do not depend on 
the relaminarization process properties and may be used 
for the arbitrary reinjection probability Pin {x) . 

3 Numerical results 

To verify the obtained theoretical predictions, we consider 
numerically the intermittent behavior of the quadratic 
map ([2|) with a stochastic force. The reinjection proce- 
dure has been fulfilled as follows: when the value of the 
variable leaves the interval — 1 < x < -1-1, the next its 
value has been taken as Xn+i = lO^^Xn- Since the depen- 
dence of the mean laminar phase length on the criticality 
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Fig. 3. The distribution of the laminar phase lengths for 
map Q. The criticality parameter value has been selected as 
e = -10"^ the noise intensity is _D = 10 ^ . The ordinate axis 
is presented in the logarithmic scale. The theoretical exponen- 
tial law (|25|l is shown by the solid line 

parameter ([^0]) has already been studied [TS], in our cal- 
culations we focus on the consideration of distribution of 
the laminar phase lengths. 

Obviously, if the intensity of noise D is equal to zero 
the intermittent behavior is observed for e > 0, whereas 
the stable fixed point takes place for e < 0. Having added 
the stochastic force we can expect that the intermittent 
behavior may be also observed in the area of the negative 
values of the criticality parameter e. 

The distribution of the laminar phase lengths is in very 
good accordance with the exponential law (j25p predicted 
by the theory of the type-II intermittency with noise (see 
Fig. [3]). Note the presence of the small region of the short 
laminar phase lengths in Fig. [3] where the deviation from 
the prescribed exponential law (1251) is observed. This re- 
gion corresponds to the transient < i < tt,- when the 



probability density px{x,t) evolves to the form (fT5)l be- 
ing close to stationary as it was discussed above. The ex- 
istence of this transient time interval does not influence 
practically on the characteristics (j20p and ((25)) of the in- 
termittent behavior in the presence of noise in the full 
agreement with the conclusions made above. So, the in- 
termittent behavior observed in the quadratic map with 
the stochastic force agrees well with the theoretical pre- 
dictions. 

4 Conclusion 

In conclusion, we have reported a type of intermittency 
behavior caused by the cooperation between the determin- 
istic mechanisms and random dynamics. The distribution 
of the laminar phase lengths being one of the important 
characteristics of the intermittent dynamics has been de- 
duced analytically. Having considered the standard model 
of type-H intermittency in the presence of noise we can 
conclude that (i) noise induces new features in the inter- 
mittent behavior of a system demonstrating type-H inter- 
mittency, with new dynamical properties being observed 
above the former value of the criticality parameter; (ii) the 
results of numerical simulations are in excellent agreement 
with the developed theory; (iii) the relaminarization pro- 
cess properties and the reinjection probability do not seem 
to play a major role for the statistical characteristic of 
type-H intermittency, and, the obtained expression of the 
mean length of the laminar phases on the criticality pa- 
rameter as well as the dependence of the mean length of 
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the laminar phases on the criticality parameter remain 
correct for different forms of the reinjection probability. 

Though the characterization of the intermittent pro- 
cess has been explicitly derived here for model system, we 
expect that the very same mechanism can be observed in 
many other relevant circumstances where the level of nat- 
ural noise is sufficient, e.g. in the physiological [^[^[^ 
or physical systems 0. 
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